Multi-body simulation of a canine hind limb: model development, experimental validation and calculation of ground reaction forces

Background: Among other causes the long-term result of hip prostheses in dogs is determined by aseptic loosening. A prevention of prosthesis complications can be achieved by an optimization of the tribological system which finally results in improved implant duration. In this context a computerized model for the calculation of hip joint loadings during different motions would be of benefit. In a first step in the development of such an inverse dynamic multi-body simulation (MBS-) model we here present the setup of a canine hind limb model applicable for the calculation of ground reaction forces. Methods: The anatomical geometries of the MBS-model have been established using computer tomography- (CT-) and magnetic resonance imaging- (MRI-) data. The CT-data were collected from the pelvis, femora, tibiae and pads of a mixed-breed adult dog. Geometric information about 22 muscles of the pelvic extremity of 4 mixed-breed adult dogs was determined using MRI. Kinematic and kinetic data obtained by motion analysis of a clinically healthy dog during a gait cycle (1 m/s) on an instrumented treadmill were used to drive the model in the multi-body simulation. Results and Discussion: As a result the vertical ground reaction forces (z-direction) calculated by the MBS-system show a maximum deviation of 1.75%BW for the left and 4.65%BW for the right hind limb from the treadmill measurements. The calculated peak ground reaction forces in z- and y-direction were found to be comparable to the treadmill measurements, whereas the curve characteristics of the forces in y-direction were not in complete alignment. Conclusion: In conclusion, it could be demonstrated that the developed MBS-model is suitable for simulating ground reaction forces of dogs during walking. In forthcoming investigations the model will be developed further for the calculation of forces and moments acting on the hip joint during different movements, which can be of help in context with the in silico development and testing of hip prostheses.

The Milky Way's Dark Matter Halo Appears To Be Lopsided

The atomic hydrogen gas (H I) disk in the outer region (beyond similar to 10 kpc from the center) of Milky Way can provide valuable information about the structure of the dark matter halo. The recent three-dimensional thickness map of the outer H I disk from the all sky 21 cm line Leiden/Argentine/Bonn survey, gives us a unique opportunity to investigate the structure of the dark matter halo of Milky Way in great detail. A striking feature of this new survey is the north-south (N-S) asymmetry in the thickness map of the atomic hydrogen gas. Assuming vertical hydrostatic equilibrium under the total potential of the Galaxy, we derive the model thickness map of the H I gas. We show that simple axisymmetric halo models, such as softened isothermal halo (producing a flat rotation curve with V-c similar to 220 km s(-1)) or any halo with density falling faster than the isothermal one, are not able to explain the observed radial variation of the gas thickness. We also show that such axisymmetric halos along with different H I velocity dispersion in the two halves, cannot explain the observed asymmetry in the thickness map. Amongst the nonaxisymmetric models, it is shown that a purely lopsided (m = 1, first harmonic) dark matter halo with reasonable H I velocity dispersion fails to explain the N-S asymmetry satisfactorily. However, we show that by superposing a second harmonic (m = 2) out of phase onto a purely lopsided halo, e. g., our best fit and more acceptable model A (with parameters epsilon(1)(h) = 0.2, epsilon(2)(h) = 0.18, and sigma(H I) = 8.5 km s(-1)) can provide an excellent fit to the observation and reproduce the N-S asymmetry naturally. The emerging picture of the asymmetric dark matter halo is supported by the. cold dark matter halos formed in the cosmological N-body simulation.

Sistema de multiprocessamento para simulação de N-corpos em FPGA

Trabalho Final de Mestrado para obtenção do grau de Mestre em Engenharia de Electrónica e Telecomunicações

IaaS clouds vs. clusters for HPC: a performance study

The increasing amount of data collected in the fields of physics and bio-informatics allows researchers to build realistic, and therefore accurate, models/simulations and gain a deeper understanding of complex systems. This analysis is often at the cost of greatly increased processing requirements. Cloud computing, which provides on demand resources, can offset increased analysis requirements. While beneficial to researchers, adaption of clouds has been slow due to network and performance uncertainties. We compare the performance of cloud computers to clusters to make clear the advantages and limitations of clouds. Focus has been put on understanding how virtualization and the underlying network effects performance of High Performance Computing (HPC) applications. Collected results indicate that performance comparable to high performance clusters is achievable on cloud computers depending on the type of application run.

Dynamische Standsicherheit von Portalstaplern

Der Artikel beschreibt die Untersuchung der dynamischen Standsicherheit von Portalstaplern. Dazu wurde ein aktueller repräsentativer Portalstapler analysiert und in ein Mehrkörper-Simulationsmodell abgebildet. Die Validierung des Modells erfolgte auf Basis von Messfahrten am realen System. Aus den durchgeführten Simulationsstudien wurden technische Verbesserungsmaßnahmen abgeleitet, welche die dynamische Standsicherheit von Portalstaplern erhöhen.

Simulationsgestützte Untersuchung der Spurtreue von Routenzügen

Für den innerbetrieblichen Transport und besonders für die Materialversorgung in der Produktion werden zunehmend Routenzüge eingesetzt, die aus einem Schleppfahrzeug und bis zu fünf Anhängern bestehen. Um gute Manövrierbarkeit und einen geringen Platzbedarf zu gewährleisten, sollten Routenzuganhänger möglichst spurtreu sein. In diesem Beitrag wird die Spurtreue exemplarisch für zwei am Markt verfügbare Anhänger für Routenzüge mit zwei ungelenkten und vier gelenkten Rädern untersucht. Hierfür wird zunächst ein Gütekriterium definiert, das die maximale Abweichung von der Spurtreue quantitativ erfasst. Außerdem werden Testszenarien vorgeschlagen, um die verschiedenen Fahrwerks- und Lenkkonzepte vergleichen zu können. Mit Hilfe eines entwickelten analytischen Modells werden die Abweichungen in der Spurtreue bei stationärer Kreisfahrt für die zwei gewählten Konzepte berechnet und dargestellt. Zusätzlich wird eine Mehrkörper-Simulation, die eine tiefere physikalische Modellierung und die Untersuchung komplexerer Fahrmanöver erlaubt, durchgeführt und mit den analytischen Ergebnissen verglichen. Es kann gezeigt werden, dass neben der Lenkkinematik weitere Parameter wie das Schräglaufverhalten der Reifen Einfluss auf die Spurtreue haben.

Domain decomposition hybrid method for numerical simulation of bluff body flows:theoretical model and application

The discrete vortex method is not capable of precisely predicting the bluff body flow separation and the fine structure of flow field in the vicinity of the body surface. In order to make a theoretical improvement over the method and to reduce the difficulty in finite-difference solution of N-S equations at high Reynolds number, in the present paper, we suggest a new numerical simulation model and a theoretical method for domain decomposition hybrid combination of finite-difference method and vortex method. Specifically, the full flow. field is decomposed into two domains. In the region of O(R) near the body surface (R is the characteristic dimension of body), we use the finite-difference method to solve the N-S equations and in the exterior domain, we take the Lagrange-Euler vortex method. The connection and coupling conditions for flow in the two domains are established. The specific numerical scheme of this theoretical model is given. As a preliminary application, some numerical simulations for flows at Re=100 and Re-1000 about a circular cylinder are made, and compared with the finite-difference solution of N-S equations for full flow field and experimental results, and the stability of the solution against the change of the interface between the two domains is examined. The results show that the method of the present paper has the advantage of finite-difference solution for N-S equations in precisely predicting the fine structure of flow field, as well as the advantage of vortex method in efficiently computing the global characteristics of the separated flow. It saves computer time and reduces the amount of computation, as compared with pure N-S equation solution. The present method can be used for numerical simulation of bluff body flow at high Reynolds number and would exhibit even greater merit in that case.

A simulation model of body contact structure in PD SOI analogue circuit

As a solution of accurate simulation of the body effect in PD SOI analogue circuit, a simulation model of distributed body contact resistance and parasitical capacitance is presented. Based on this model, we have designed and simulated a sense amplifier that applied to V a 0.8um PD SOI 64K SRAM.

Molecular dynamics simulation of fractures using an N-body potential

Molecular dynamics has been employed to model the fracture of a twodimensional triangular atomic lattice. The N-body Sutton-Chen potential developed for fcc metals and its extended version (Rafii-Tabar and Sutton) for fcc random binary alloys were used for the interatomic interactions. It is shown that at low temperatures cleavage fractures can occur in both an elemental metal and an alloy. At elevated temperatures the nucleation of dislocations is shown to cause a brittle-to-ductile transition. For the brittle crack propagation in the elemental metal, crack propagation speeds have been computed for different stress rates, and a crack instability found to exist as the speed reaches a critical value of about 32% of the Rayleigh wave speed. For the random alloy, we find that the dislocation movement can be affected by the distorted lattice.

Molecular dynamics simulation of fractures using an N-body potential

Molecular dynamics has been employed to model the fracture of a two dimensional triangular atomic lattice. The N-body Sutton-Chen potential developed for fcc metals and its extended version (Rafii-Tabar and Sutton) for fcc random binary alloys were used for the interatomic interactions. It is shown that at low temperatures cleavage fractures can occur in both an elemental metal and an alloy. At elevated temperatures the nucleation of dislocations is shown to cause a brittle-to-ductile transition. For the brittle crack propagation in the elemental metal, crack propagation speeds have been computed for different stress rates, and a crack instability found to exist as the speed reaches a critical value of about 32% of the Rayleigh wave speed. For the random alloy, we find that the dislocation movement can be affected by the distorted lattice.

Measurements, modeling and simulation of the off-body radio channel for the implementation of bodyworn antenna diversity at 868 MHz

This paper presents a systematic measurement campaign of diversity reception techniques for use in multiple-antenna wearable systems operating at 868 MHz. The experiments were performed using six time-synchronized bodyworn receivers and considered mobile off-body communications in an anechoic chamber, open office area and a hallway. The cross-correlation coefficient between the signal fading measured by bodyworn receivers was dependent upon the local environment and typically below 0.7. All received signal envelopes were combined in post-processing to study the potential benefits of implementing receiver diversity based upon selection combination, equal-gain and maximal-ratio combining. It is shown that, in an open office area, the 5.7 dB diversity gain obtained using a dual-branch bodyworn maximal-ratio diversity system may be further improved to 11.1 dB if a six-branch system was used. First-and second-order theoretical equations for diversity reception techniques operating in Nakagami fading conditions were used to model the postdetection combined envelopes. Maximum likelihood estimates of the Nakagami-parameter suggest that the fading conditions encountered in this study were generally less severe than Rayleigh. The paper also describes an algorithm that may be used to simulate the measured output of an M-branch diversity combiner operating in independent and identically-distributed Nakagami fading environments.

Modal Representation of the Resonant Body within a Finite Difference Framework for Simulation of String Instruments

This paper investigates numerical simulation of a string coupled
transversely to a resonant body. Starting from a complete nite
difference formulation, a second model is derived in which the
body is represented in modal form. The main advantage of this hybrid form is that the body model is scalable, i.e. the computational
complexity can be adjusted to the available processing power. Numerical results are calculated and discussed for simplied models
in the form of string-string coupling and string-plate coupling.

Many-body dissipative particle dynamics simulation of wetting phenomena

With the development of the simulation of particle dynamics, the traditional dissipative particle dynamics (DPD) method can not satisfy the needs of research in static or dynamic wetting phenomena. However, the Many-body DPD approach extends the ability of the traditional method to simulate the interface between solid and liquid or some other situation. In this paper, we propose a Many-body DPD program to simulate the solid-liquid interface and get satisfactory results.

A full body musculoskeletal model based on flexible multibody simulation approach utilised in bone strain analysis during human locomotion

Load-induced strains applied to bone can stimulate its development and adaptation. In order to quantify the incident strains within the skeleton, in vivo implementation of strain gauges on the surfaces of bone is typically used. However, in vivo strain measurements require invasive methodology that is challenging and limited to certain regions of superficial bones only such as the anterior surface of the tibia. Based on our previous study [Al Nazer et al. (2008) J Biomech. 41:1036–1043], an alternative numerical approach to analyse in vivo strains based on the flexible multibody simulation approach was proposed. The purpose of this study was to extend the idea of using the flexible multibody approach in the analysis of bone strains during physical activity through integrating the magnetic resonance imaging (MRI) technique within the framework. In order to investigate the reliability and validity of the proposed approach, a three-dimensional full body musculoskeletal model with a flexible tibia was used as a demonstration example. The model was used in a forward dynamics simulation in order to predict the tibial strains during walking on a level exercise. The flexible tibial model was developed using the actual geometry of human tibia, which was obtained from three-dimensional reconstruction of MRI. Motion capture data obtained from walking at constant velocity were used to drive the model during the inverse dynamics simulation in order to teach the muscles to reproduce the motion in the forward dynamics simulation. Based on the agreement between the literature-based in vivo strain measurements and the simulated strain results, it can be concluded that the flexible multibody approach enables reasonable predictions of bone strain in response to dynamic loading. The information obtained from the present approach can be useful in clinical applications including devising exercises to prevent bone fragility or to accelerate fracture healing.